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Abstract. Weakly interacting massive particles (WIMPs) are arguably the most 
natural DM candidates from a particle physics point of view. After their number 
density has frozen out in the early universe, thereby setting their relic density today, 
WIMPs are still kept very close to thermal equilibrium by scattering events with 
standard model particles. The complete decoupling from the thermal bath happens as 
late as at temperatures of around ~ 1 — 10 MeV and provides an important cosmological 
scale that can directly be translated into a small-scale cutoff of the spectrum of 
matter density fluctuations. We present here a full analytic treatment of the kinetic 
decoupling process from first principles. This allows an exact determination of the 
decoupling scale, for arbitrary WIMP candidates and any scattering amplitude. As 
an application, we consider the situation of the lightest supersymmetric particle as 
well as the lightest Kaluza-Klein particle that arises in theories with universal extra 
dimensions; furthermore, we show that our formalism can also easily be applied to, e.g., 
the evolution of the non-relativistic electrons into the recombination regime. Finally, 
we comment on the impacts for the smallest gravitationally bound structures and the 
prospects for the indirect detection of dark matter. 
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1. Introduction 

Based on the simplest inflationary scenario, a flat cosmological model with dark energy 
and cold dark matter as dominant contingents in the cosmic energy budget fits the 
spectrum of observed temperature anisotropies imprinted in the cosmic microwave 
background radiation (CMB) pQ. Assuming that the primordial fluctuations are 
adiabatic with a power law spectrum, the temperature anisotropy data alone indicates a 
significant amount of non-baryonic cold dark matter (CDM); analysis of the anisotropies 
finds that the relative matter density u m = Q m h 2 = 0. 127^0 013 is roughly six times larger 
than the relative baryon density. This is consistent with a diversity of other cosmological 
observations, e.g. the analysis of galaxy clustering as observed by SDSS [2] and 2dFGRS 
[3], improvements in the lensing data [I], as well as improvements in small-scale CMB 
measurements, e.g. [5]. 

Arguably the best studied class of CDM particle candidates are weakly interacting 
massive particles, or WIMPs (for reviews on particle DM, see [6j). Equipped with 
masses and typical interaction strengths set by the electroweak scale, they naturally 
occur in many extensions to the standard model of particle physics (SM), e.g. in its 
super symmetric versions or in models with space-like extra dimensions. Thermally 
produced in the early universe, and provided their stability, WIMPs automatically 
acquire a relic abundance that is of the right order of magnitude to account for the 
observed matter density. Another well-studied CDM particle candidate is the axion. 
The essential difference between WIMPs and axions is that the former have been kept 
in local thermal equilibrium prior to some epoch in the cosmological evolution while the 
latter have not. This is because by the time the axion potential developed a minimum at 
around the QCD phase transition, and the axions thus became massive, any interactions 
were already heavily suppressed by the Peccei-Quinn scale. In this article, we will focus 
solely on WIMPs and often use the abbreviation CDM as a synonym for these. 

It has long been a widely held belief that the process of structure formation is 
completely independent of the CDM particle nature; CDM density perturbations, in 
particular, should only be effected by their gravitational interaction and bear no visible 
imprint of particle interactions. In this case, however, cosmological perturbation theory 
predicts a monotonously growing support for structure formation on ever smaller scales. 
More precisely, the spectrum of matter density perturbations grows logarithmically 
with the comoving wavenumber that characterizes the perturbation mode [3, E] • This 
behaviour has recently been confirmed yet another time, down to the resolution 
of a numerical simulation with unprecedented accuracy [14] . As a consequence, a 
regularization scheme would be required to render linear perturbation theory well- 
defined, e.g. a cut-off in the wavenumber that corresponds to some resolution limit. 
From cosmological perturbation theory we know that gravitational interactions alone 
cannot give rise to a fundamental cut-off scale, so it has to be introduced by hand. On 
the other hand, it is of course rather problematic to impose such an ad hoc small-scale 
cut-off since, in hierarchical structure formation, small scale perturbations enter the 
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nonlinear regime first and structures on larger scales form only later via mergers and 
accretion. 

In Ref. [9] it was shown that CDM interactions with the radiation component give 
rise to a natural cut-off that can be calculated directly from the interaction theory 
underlying a specific dark matter candidate. The authors of pU] confirmed the results 
of [9] . The (physical) cut-off is set by viscous processes around the kinetic decoupling of 
CDM from radiation. Close to kinetic decoupling, CDM density perturbations receive 
in-medium modifications that give rise to a dispersion relation with a non-vanishing 
imaginary contribution and thus to absorption. Subsequently, free streaming of CDM 
suppresses further the power of structure formation on small scales. In that way, 
collisional and collision-less damping processes give rise to a fundamental cut-off in 
the spectrum of linear CDM perturbations [7J, [8] and thus to a characteristic feature 
that clearly depends on the CDM micro-physics. It was later shown numerically [TT] 
that acoustic oscillations have to be taken into account at kinetic decoupling, however, 
the order of magnitude of the cut-off scale is correctly given in [71 [8]. Most recently, the 
author of [12] succeeded in including acoustic oscillations before and during kinetic 
decoupling, the viscous coupling during kinetic decoupling and free streaming in a 
consistent kinetic treatment. Close to this physical cut-off scale, the power spectrum 
develops a maximum. Hence, this scale characterizes those density perturbations that 
typically enter the non-linear regime first and therefore give rise to the first purely 
gravitationally bound structures. 

Provided that the dark matter density distribution today retains traces of 
these typical first halos, the cut-off scale described above may thus have important 
observational consequences. While it is certainly true that the present CDM density 
profile is independent of the particle nature of CDM on super-kpc scales, this is 
highly questionable on sub-pc scales. Let us therefore stress that WIMP detection 
experiments actually probe the dark matter density profile on sub-mpc scales [13J: for 
a given WIMP candidate, the prospects for a direct detection (via elastic scattering 
processes with the detector nuclei) depend on the effective reaction volume which is 
given by the detector size and the path the detector travels during data taking. Indirect 
detection experiments, on the other hand, aim at measuring flux densities from dark 
matter annihilation products, which depend quadratically on the dark matter density 
contrast in the observed region. Substructures thus enhance the fluxes and therefore the 
prospects for indirect detection, probing the dark matter density profile on sub-galactic 
scales. 

Let us stress that it is, in fact, the temperature at kinetic decoupling that provides 
the interface, at least in principal, between the interaction theory underlying a specific 
CDM candidate and CDM substructures that are probed by direct and indirect dark 
matter searches. Motivated by this observation, we present in this article a complete 
and analytic description of the kinetic decoupling of WIMPs from the local heat bath, 
without assuming a specific WIMP candidate. In Section [2j we introduce important 
time scales that characterize the non-equilibrium state of WIMPs and discuss on a very 
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qualitative level both chemical and kinetic decoupling of WIMPs, before we turn to 
the presentation of a precise treatment of the kinetic decoupling process. We introduce 
a generalized WIMP temperature (which characterizes the WIMP distribution even in 
the kinetic stage and agrees with the heat bath temperature before kinetic decoupling) 
and derive a master equation for the the evolution of this temperature. Section [3] is 
devoted to various examples: as a first and familiar example we quantify, for illustration, 
the deviation of the electron temperature from the photon temperature during the 
recombination regime. More relevant for our original purposes is the next example, 
the kinetic decoupling of bino-like CDM from SM fermions. As a last example we 
consider Kaluza-Klein dark matter. Section H] provides a summary of the processes that 
constitute the small-scale cut-off in the linear matter power spectrum; here, we also 
shortly review what is known about the total mass fraction in and survival probability of 
the smallest gravitationally bound structures that correspond to this scale. In section [5] 
we discuss prospects for indirect dark matter detection experiments in light of small-scale 
structure formation. Section El finally, concludes. All technical details and derivations 
are presented in the three appendices: In Appendix A , we solve the Boltzmann equation 
for the collision integral, as determined in Appendix B , that describes the scattering 
processes of relevance to our discussion; in Appendix C , we present the LKP scattering 
amplitude and use this opportunity to demonstrate the applicability of our formalism 
even in the case of s-channel resonances of the scattering amplitude. 



2. Thermal relics and their decoupling from the heat bath 

In this section, we introduce various time scales that allow to characterize different 
stages in the evolution of WIMPs. While a thorough understanding of the various 
characteristic time scales is not required for the later computation of the kinetic 
decoupling temperature, the presentation mainly aims at clarifying the importance of the 
relaxation time for WIMPs close to local thermal equilibrium. We then discuss chemical 
and kinetic decoupling of WIMPs from the local heat bath on a very qualitative level, 
before proceeding to a precise description of kinetic decoupling - which is the main 
concern of this article. 



2.1. Characteristic time scales 

The most detailed statistical description of WIMPs is given by the iV-particle phase 
space density, with N denoting their number in a comoving Volume V. Around kinetic 
decoupling, the interaction radius r- m t of WIMPs (defined as the maximal distance 
between a WIMP and SM - or, rather, heat bath - particles over which interactions can 
occur) is much smaller than the mean distance between the particles, i.e. Nr? t /V <C 1. 
This simple inequality allows the introduction of different time scales that are important 
for the classical kinetic description of dark matter. The first time scale is called 
collision time and defined by r c = ri n t/(v), with (v) denoting the mean dark matter 
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velocity. A second time scale is provided by the time between two subsequent collisions, 
Tf = If I (v) , where If is the mean free path of dark matter particles. The third and 
for our considerations most important time scale is the relaxation time r r . It marks 
the time during which local thermal equilibrium is established in a comoving volume 
that contains many dark matter particles, but that is still small compared to V. The 
last time scale we want to consider here, r eq , is given by the time during which global 
equilibrium in V is established. 

From elementary kinetic theory it follows that the mean free path can be 
approximated by lf/V 1/3 ~ (1/N)(V 1/3 /r int ) 2 , hence r c < r f . Since CDM particles with 
weak interactions need many collisions in order to maintain local thermal equilibrium, 
we also have Tf <C r r . Finally, one evidently has r r <C r cq . To summarize, we get the 
following chain of inequalities for the characteristic time scales introduced above: 

r c <T f <r r <r cq . (1) 

For a description of processes on time scales 5t <C r c , the full iV-particle phase space 
density is required, and a reduced description in terms of some N < iV-particle phase 
space distribution is impossible. A reduced kinetic description is, however, available 
for processes with r c <C 5t <C r r . The WIMPs can then adequately be described by a 
one-particle distribution that determines the probability distribution for the locations 
and the momenta of the CDM particles. For processes on time scales r r <C 5t <C r eq , 
the system is in local thermal equilibrium and further simplifications apply. During 
this stage, the WIMPs are fully described by the local number, momentum and kinetic 
energy density, which are given as the zeroth, first and second moments, respectively, 
of the one-particle distribution with respect to the WIMP momentum. This is the 
underlying reason that allows to describe such a system as a perfect or almost perfect 
fluid. 

2. 2. Scales of decoupling 

Initially, the WIMP system is (almost) in local chemical and thermal equilibrium, with 
chemical equilibrium provided by the detailed balance between WIMP annihilation and 
production processes, and local thermal equilibrium maintained by elastic scattering 
processes with SM particles in the heat bath. Typically, chemical and thermal 
equilibrium are protected by the aforementioned processes for temperatures above the 
WIMP particle mass, T > M x . 

As the Universe expands, the temperature decreases and the WIMP momenta 
become redshifted. At T ~ M x detailed balance is destroyed: the lighter species, 
dominating the energy budget, can no longer produce WIMPs, while the reverse channel 
is still open. As a consequence, between the end of detailed balance at T ~ M x and 
chemical decoupling at T = T C( j, the total (comoving) number of WIMPs decreases. 
Chemical decoupling happens roughly when the WIMP annihilation rate drops below the 
Hubble rate H, and, as a consequence, annihilation processes cease to proceed and the 
relic dark matter abundance is fixed. This is expected to happen since the relevant target 
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density in the annihilation processes is provided by the WIMP number density, which 
is Boltzmann suppressed for temperatures T < M x . For instance, in the case of WIMPs 
the chemical decoupling scale is typically given by T cd ~ M x /25. Hence, the total 
number of WIMPs is conserved for T > M x - dynamically protected by detailed balance, 
decreases during M x /25 < T < M x , and again protected for T < M x /25 - basically by 
the volume expansion]]] Note that the dependence of the chemical decoupling scale T cd 
on the WIMP annihilation rate (and thus on the underlying interaction theory) can be 
completely eliminated by using the value for the relative WIMP density as inferred e.g. 
from the temperature anisotropy data. 

After chemical decoupling elastic scattering processes still keep WIMPs in local 
thermal equilibrium. It is very important to realize that even though dark matter 
is now a chemically distinct particle species, it is still a constituent of the local heat 
bath. The relevant target density for elastic scattering processes is provided by the 
number density of relativistic SM particles and thus only decreasing as T 3 . Eventually, 
also the elastic scattering rate r e i drops below the Hubble rate and elastic scattering 
processes cease. However, even before last scattering, local thermal equilibrium cannot 
be maintained since the relaxation processes within the WIMP system are weakened by 
the Hubble expansion. Consequently, the relevant time scales are the relaxation time r r 
(intrinsic) and the Hubble time H~ l (extrinsic). 

Before we will present a way to analytically determine the decoupling scale in the 
next Section, let us now roughly estimate the relaxation time for WIMPs as follows 
(see also [9]): our system consists of a weakly interacting particle species with mass 
M x and a local heat bath characterized by a temperature T -C M x , so M x is the 
dominant energy scale. The average (three-) momentum transfer Ap from a light (heat 
bath) particle to the WIMP, during a single elastic scattering process, is thus roughly 
given by its energy, i.e. Ap ~ T. We assume now that Ap is constant during the 
relaxation process. After N co \\ collisions, inducing a random walk in momentum space, 
the WIMP momentum will thus have changed by Ap (iV co n) 1//2 . In order to maintain local 
thermal equilibrium, this should be comparable to the typical momentum for thermally 
distributed WIMPs, which is given by (M x T) l l 2 . The number of collisions required is 
thus iV co ii ~ M x /T ^> 1 and the relaxation time becomes r r ~ iV co ii/r e i ~ (m/T) l/<7 sca i e . 
The scale for the total elastic scattering process can furthermore be assumed to be 
approximated by cr scale = (Gpm^M^/m^, with Gp denoting Fermi's constant, m w and 
mz are the masses of the charged and neutral weak gauge bosons, respectively. Kinetic 
decoupling occurs approximately at the temperature T kd , for which r r (T kd ) = if _1 (T kd ). 
With these simplifications, we thus expect T kd ~ M x ^ A mz/M^( i {G-F , m^ ! ) 1 / 2 ~ (9(1 — 10) 
MeV. Let us stress again that this is only a rough estimate for the decoupling scale; 
below, we will determine the precise answer by solving the appropriate kinetic theory. 

% Strictly speaking, the relic abundance is protected only until density perturbations generate density 
contrasts that again allow for WIMP annihilations. 
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2. 3. A precise description of the kinetic decoupling process 

Around kinetic decoupling WIMPs pass from the hydrodynamical to the kinetic stage. 
Their appropriate description changes accordingly: while the hydrodynamical stage 
describes the evolution of fluid variables that are given as expectation values with 
respect to the one-particle phase space distribution, the kinetic stage is concerned with 
the evolution of the one-particle phase space distribution /(p) itself. 

A more precise determination of the decoupling temperature requires the 
introduction of a temperature parameter T x for dark matter, which we choose to define 
as the second moment of WIMP phase-space distribution function /: 

r H(i)' with <°< p » s U °™ • < 2 » 

Here, n x denotes the WIMP number density and O is some microscopic observable. 
The parameter T x that we introduced above agrees with the heat bath temperature 
for T > Tkd and can be viewed as a generalized WIMP temperature for T < its 
difference to T thus characterizes the deviation of WIMPs from the equilibrium state. 

The dynamics of the one-particle phase space distribution, governed by the 
Boltzmann equation, is nontrivial during the kinetic decoupling process (kinetic stage). 
In the case of dark matter weakly coupled to radiation the main complication arises 
from the coupling of slow processes (non-relativistic propagation of dark matter) and 
fast processes (relativistic propagation of radiation) in the collision integral, which is 
enforced via the explicit energy momentum conservation constraining the phase space 
available for elastic scattering processes. 

Around kinetic decoupling, the WIMP mass M x is by far the largest energy scale 
involved. In particular, the kinetic WIMP energy, as well as the average momentum 
transfer during the elastic scattering processes, is tiny compared to M x . This allows to 
factorize slow and fast processes in the collision integral by utilizing a Born-Oppenheimer 
type of approximation and expanding the collision integral in p 2 /M x 2 . Instead of solving 
the Boltzmann equation for /, we can then solve its second moment and derive, in 
Appendix A, the following process equation for T X (T): 



dT 



2+a[ k 



T x ( T 



M x 



(3) 



Here, n is determined by the low-energy behaviour of the scattering amplitude at zero 
momentum transfer, |A^|^ =0 cx (uj/M x ) n , where u is the energy of the SM scattering 
partner; a is a numerical constant that is defined in (A. 23). The uniquely defined 
solution to Eq. (j3J) with the correct asymptotic behavior in the high-temperature limit, 
T x T ^°° T, is given by: 

f rV(«+2) ^ 

T * = T { 1 - "STT^l r H» + 2) • 4 L ' (4) 



n + 2 V M- 



In the low-temperature limit, T — > 0, we have T x ~ T 2 /M x . Since, as we shall see, 
kinetic decoupling occurs on a rather short time scale, a suitable choice for a definition 
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of the kinetic decoupling temperature is obtained by matching the two above-mentioned 
asymptotic solutions to Eq. (j3J), resulting in: 

l/(n+2) 



T kd ( T 



M x \M X T X J T ^ \\n + 2 



r 



n + l 
n + 2 



(5) 



For a more detailed discussion of the process equation ([3]), we refer to Appendix A 



and the following Section, where we consider some examples of particular interest and 
determine the exact values of the decoupling temperature in these cases. However, from 
( 1A.23I) . let us already now remark that we generically expect it to be of the order of 



(a-*M x /M Pi y<^> , (6) 



where a is the effective coupling constant that appears in the scattering amplitude. For 
example, we have a = a em , n = in the case of recombination (Thompson scattering) 
and a = ay, n = 2 for the thermal decoupling of Binos or Kaluza-Klein photons. From 
Eq. ([6]) we can see that, indeed, the WIMP micro-physics enters in a crucial way in the 
determination of the kinetic decoupling scale and thus - see Section H]- in the numerical 
value for the cutoff in the spectrum of CDM density perturbations. 



3. Examples 

In this section, we apply the formalism developed in Appendix A and B, as reviewed in 
the preceding Section, for illustrative purposes first to the decoupling of electrons from 
radiation during the recombination era. In a second step, we turn closer to the main 
focus of this article and discuss the kinetic decoupling of typical WIMP DM candidates, 
i.e. bino-like and Kaluza-Klein dark matter, in particular. 



3.1. Recombination 

After the (chemical) decoupling of electrons and positrons at a temperature of around 
T ~ 1 MeV, and the subsequent annihilation of e + e~ pairs, the photons were 
still interacting with protons and the remaining electrons, mainly through Compton 
scattering with the latter. Such a plasma of interacting relativistic and non-relativistic 
particles cannot, strictly speaking, be kept in thermal equilibrium in an expanding 
universe. On the other hand, CMB observations show a black-body spectrum to an 
accuracy of about 10~ 4 [15J. As a first application of our formalism, we will now quantify 
the expected deviation from a thermal distribution at the time when the universe became 
transparent and the CMB photons that we observe today were released. For a similar 
treatment of the transition into the recombination regime, see, e.g. [16J. 

At zero momentum transfer and for small fermion energies, the amplitude for 
Thompson scattering is given by 

\M\ 2 t=0 = 6^ 2 a 2 em . (7) 
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From our general solution (111) , the temperature of the non-relativistic electrons therefore 
evolves as 

T e = T 1- exp[z]r [-1/2, z}\ (8) 
= Tv^riexp[z] {l-Erf [Vf|}^B/_T.\a , (9) 

Z -2\m e J 



where 



1X231 / 10 ino 2 2 Ar _Mpi 



(2T).. 3.36 128^,^^ = 1-8 x 10» . (10) 

We thus find that when entering the recombination regime at around T ~ 0.3 eV, 
the electron temperature deviates only by (T — T e )/T = 1.6 x 10~ 7 from the photon 
temperature. The electrons, in other words, do not in any significant way affect the 
photon temperature and the photons thus follow an equilibrium distribution into the 
recombination regime, explaining the observed black-body nature of the CMB. This 
argument can be made even more quantitative by considering the Boltzmann equation 
for the photons which, in this regime, is known as the Kompaneets equation [T7j : 

(d t - Hud u ) g(u, t) = *^ n ^d u [u*{T e d u g(u, t) + [l + g(u>, t)] g(u, t)}] , (11) 

where g and u are the photon distribution function and energy, respectively. One may 
now define the photon temperature by the relation 

g(l + g) = -Td w g, (12) 

which is automatically satisfied for a thermal distribution. With this definition, one 
realizes that for T = T e the collision term on the right hand side of fTTTl) vanishes 
identically and the photons therefore just follow an equilibrium distribution. 



3.2. Neutralino dark matter 

One of the most well-motivated and extensively studied DM candidates is the neutralino. 
In most supersymmetric extensions to the SM, it appears as the lightest supersymmetric 
particle (LSP) and is given by a linear combination of the super-partners of the gauge 
and Higgs fields, 

X = X°i = N X1 B + N 12 W 3 + iV^ + iV 14 ^2° • (13) 

Here, in order to illustrate our formalism and to describe the kinetic decoupling process 
for a typical WIMP dark matter candidate, we consider the simple case of a Bino-like 
neutralino, x ~ B- From the WMAP relic density constraint, Binos have a mass of a 
few hundred GeV or less. 

To lowest order in the fermion energy u, the scattering amplitude of Binos with 
SM fermions is given by [9]: 

\M\U = 32*4 M + n 4 ) f - 2 • (14) 
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Figure 1. The LSP temperature T x as a function of the background temperature T in 
the early universe is shown as a solid line, for the case of a Bino with Ms = 100 GeV 
and Mi = 150 GeV. The rather abrupt transition to a regime with T x oc T 2 , where 
the LSP is completely decoupled from the thermal bath, makes it straightforward to 
identify the temperature Tkd of thermal decoupling as given in Eq. (| 18[) : here, it is 
indicated by a dashed line. For comparison, the (green and blue) dotted lines show 
the naive estimates for the evolution of T x , assuming a sudden decoupling from the 
thermal bath at the time of last scattering or the relaxation time, respectively (see also 



From (j4]), we thus find that the neutralino temperature evolves with the background 
temperature as 

T x /T=\l-—exp[z]r[-l/4 c ,z]\ , (15) 



Dl' 



where 



JX23J 31 /57r 3 M P1 4 ^ 4 4 | | 





This corresponds to a decoupling temperature of 

Tkd i, 5x K_ 1 V / V 77 i^V /1 MeV. (18) 



M% I V 100 GeV 

When computing the sum over all contributions, we have assumed an identical mass for 
all relevant sleptons L (i.e. (anti-)selectrons and (anti-)sneutrinos). 

In Fig. [1} we show the evolution of the CDM temperature for a typical Bino DM 
candidate during the process of kinetic decoupling. For reference, we also indicate the 
naive estimates for the decoupling temperature that were reviewed in Section [2j At 
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temperatures only slightly greater than T^, the LSP is still very tightly bound to the 
background temperature T (for reference, see also Fig. lAlj) . We therefore expect our 
result to be highly independent of the details of the QCD phase transition at around 
Tqcd ~ 150 MeV. A second important observation is the rather sharp turnover from 
T x pa T to a regime with T x oc T 2 , where the LSP has completely decoupled from the 
thermal bath and simply cools down due to the Hubble expansion. Finally, to take 
up our previous discussion on the various scales that are involved in the decoupling 
process, we can clearly see that the relevant quantity is not the time of last scattering, 
which happens long after the LSP has already decoupled. The relaxation time r re i, on 
the other hand, gives as expected a much better estimate for the time of decoupling. 
The reason that it underestimates the actual decoupling time is to be found in the fact 
that close to decoupling only small deviations from a thermal distribution need to be 
restored - while the derivation of the relaxation time r r implicitly assumed that thermal 
equilibrium had to be restored from a highly non-thermal distribution. 

Obviously, our formalism can in principle easily be extended beyond the simple 
Bino case that we have presented here and applied to an arbitrary neutralino candidate. 
In particular, it would be interesting to perform a scan over, e.g., the MSSM parameter 
space to better understand the expected spread of the decoupling temperature. We note 
that a first step into this direction has been performed in Ref. [18], based on the estimate 
that the relaxation time as defined before is a good approximation for the decoupling 
time; in the case of the Bino, this would correspond to a kinetic decoupling temperature 
of [9] 

(m \ 1/2 / Mb \ 5/4 

^^Mf- 1 ) (lOOtfcv) ^ (19) 

which we have, for comparison, also included in Fig. [IJ We would like to stress that it 
would be warranting to re-evaluate the results of Ref. [18J in light of the method that we 
have presented here; in particular, since we are now able to actually follow the details 
of the kinetic decoupling process, the underlying uncertainty in the determination of 
the cutoff-scale is greatly reduced. Also, the method can be applied to regions of the 
parameter space where a simple scaling of the cross section with energy, like a ~ const. 
or a oc Ef, is not valid anymore (see Appendix C). 



3.3. Kaluza-Klein dark matter 

Models with universal extra dimensions (UED) [19], where all standard model (SM) 
fields are allowed to propagate in a higher-dimensional bulk, have received a great deal 
of attention since it was realized that they naturally give rise to a new class of dark 
matter candidates [20l El]: the higher- dimensional extra degrees of freedom appear 
in the low-energy effective 4D theory as towers of new, heavy states, the lightest of 
which - similar to the case of i?-parity in supersymmetry - is stable due to an internal 
Z2 symmetry ( "KK-parity" ) that appears as a remnant of the higher-dimensional 
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translational invariance; thermally produced in the early universe, the lightest Kaluza- 
Klein particle (LKP) acquires the right WMAP relic density for a compactification scale 
of about R- 1 ~ 1 TeV [221 E31 E3j. 

At tree-level, all Kaluza-Klein modes are essentially degenerate in mass, M ~ -R" 1 , 
so radiative corrections become important for the determination of the actual mass 
hierarchy between these states and, in particular, for the determination of the LKP. 
Unfortunately, the radiative mass shifts depend in principle on the (unknown) physics 
at the cutoff scale A and should in that sense be treated as free - though small - 
parameters. However, one usually follows the (self-consistent) assumption that the 
corresponding contributions at that scale are negligible [20]. In this minimal UED 
(mUED) model, all masses are determined by only two parameters, R and A, and the 
LKP turns out to be the the first KK excitation of the weak hypercharge gauge 
boson. 



In Appendix C we calculate the scattering amplitude for the LKP with SM 
fermions; to lowest order in the fermion energy lu, it is given by : 




'B(l) / \ fW BW 




(21) 



where in the last step we used the fact that the fermion singlet KK states receive much 
smaller radiative corrections to their masses than the KK doublet states and introduced 

5 S = (M Cl ) - M BW )/M BW < 1 . (22) 

We note that the amplitude takes the same form as for a Bino LSP. The evolution of 
the LKP temperature is thus described by Eq. (1151) . but now with 

31 M m 4 ^ \ A M\ m \ I M|„, V\,„, 

«2xl0»^(^i)". (24) 
The corresponding decoupling temperature is therefore given by 

T kd l3xl0 2 ^ 2 (^) 5/4 MeV. (25) 

The sum that appears above runs over all SM fermions that contribute to the scattering 
process; for temperatures around the kinetic decoupling, these are only electrons and 
the three neutrino species, as well as the corresponding antiparticles (the contribution 
from all other particles is strongly Boltzmann suppressed). The mass-splitting between 
the LKP and the corresponding KK states is then given by 5 ~ 10~ 2 . In Figure [2], we 
take the exact mUED values for 6 and show T kd as a function of R and A, using the full 
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l/R M 

Figure 2. The left part of this figure shows the kinetic decoupling temperature Tkd 
for the LKP in the mUED model. From top to bottom, the cutoff scale is set to 
AR = 40, 30, 20 and the grey band shows the region consistent with the 2a WMAP 
relic density constraint for a Higgs mass mh < 150 GeV (for higher Higgs masses, the 
grey region broadens and shifts to the right, up to < 1.3 TeV [24]). In the right 
panel, we leave the mUED setup and treat the mass splitting between the and the 
KK leptons as a single, free parameter 5; from bottom to top, the curves correspond 
to S = 0.01, 0.05, 0.2. 



expression (1231 . We find that the LKP decoupling temperature falls into a very narrow 
range, 

27MeV<T kd <34MeV, (26) 

depending on the parameters of the model. The low value of T kd , despite a TeV scale 
LKP mass, justifies the above statement about the number of relativistic SM degrees 
of freedom close to decoupling. In particular, we would like to stress that the kinetic 
decoupling of the LKP - just as for the Bino - happens sufficiently late so as not to be 
influenced by the details of the QCD phase transition at around Tqcd ~ 150 MeV. 

As a next step, we relax the assumptions of the mUED setup and treat the LKP 
mass and the mass splitting as independent parameters that may be varied about the 
mUED values. As can be seen from the right panel of Fig. [2], one may in this way 
encounter higher decoupling temperatures - but without invoking extreme fine-tuning, 
one will still stay below T kd < 100 MeV. 

4. From kinetic decoupling to damping scales in the power spectrum 

After having provided a detailed description of the kinetic decoupling process in the 
preceding part of this article, we would now like to briefly discuss in the next two 
Sections, on a rather qualitative level, the further evolution of the WIMP population in 
the expanding universe. 

To start with, and in order to get a qualitative intuition for the physical processes 
leading to kinetic damping of CDM perturbations, let us use the hydrodynamical 
description for discussing the evolution of these perturbations. For temperatures 
T > T kd , i.e. before kinetic decoupling, small-scale perturbations in the CDM fluid are 
damped due to the tight coupling between the CDM and radiation fluid. Note that this 
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includes temperatures T < M x , when the CDM fluid is already non-relativistic. Around 
kinetic decoupling, T ~ T k d, the CDM fluid starts to deviate from an ideal fluid, i.e. the 
coupling to the radiation component is not strong enough to maintain local thermal 
equilibrium for the CDM any longer. However, the residual coupling between these 
two fluids generates some friction between them. The main viscous processes causing 
this friction are bulk- and shear- viscosity, while heat conduction is sub-dominant. 
Hence, perturbations of the CDM energy momentum tensor behave like sound waves 
propagating through an opaque medium (the heat bath) [25]. Their dispersion relation 
receives in-medium modifications, resulting in an imaginary frequency that is generated 
by the viscous coupling between the CDM perturbations and the heat bath. 

The strength of the viscous couplings can not be calculated within the 
hydrodynamical framework. The most convenient strategy for calculating the viscosity 
coefficients is to calculate the entropy generation in the CDM-radiation mixture 
within kinetic theory and equate it to the corresponding expression derived in the 
hydrodynamical framework. This allows to simply read off the viscosity coefficients 
by comparing the relevant tensors in both approaches that describe, respectively, the 
same type of entropy generation. 

The effect of the viscous coupling is to damp perturbations in the CDM energy 
momentum tensor; in fact, they generate an exponential cut-off in the perturbation 
spectrum at kinetic decoupling. This spectrum then serves as initial condition for the 
subsequent evolution at temperatures T < T kd - the so-called free streaming regime, 
where a further damping in the perturbation spectrum takes place. The free streaming 
of CDM particles leaves a finger-print similar to collisional damping in the spectrum, 
i.e. a characteristic exponential damping term. Free streaming is therefore often also 
referred to as collision-less damping. For WIMPs, the characteristic damping mass is 
set by 

1 + In (T kd /30 MeV) /19.2 



M fe 
10- 6 M m 



(27) 



(M x /100 GeV) 1/2 (T kd /30 MeV) 1/2 

The characteristic damping scale for free streaming is time dependent, however, it 
becomes constant to a very good approximation sufficiently late after matter radiation 
equality. 

The authors of [TJ [5] calculated the envelope of the CDM perturbation spectrum. 
In [11] it was later shown numerically that acoustic oscillations have to be taken into 
account at kinetic decoupling as they give rise to even stronger damping; the order 
of magnitude of the cut-off scale, however, is correctly given in [8]. Most recently, 
the author of [12] succeeded in including the acoustic oscillations before and during 
kinetic decoupling, the viscous coupling during kinetic decoupling and free streaming in 
a consistent kinetic treatment. All approaches agree in that the dominant suppression of 
small-scale power is given by collision-less damping during free streaming, with kinetic 
decoupling providing consistent initial conditions. 

During radiation domination and on sub-Hubble scales, the transfer function for 
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matter density perturbations grows logarithmically, while after matter radiation equality 
matter density perturbations increase linearly (for a vanishing baryon fraction) with the 
scale factor on sub-Hubble scales. In both eras, during radiation and during matter 
domination, the spectrum of perturbations depends logarithmically on the comoving 
wavenumber that characterizes the corresponding mode. Hence, modes characterized 
by wavenumbers close to the cut-off in the power spectrum typically enter the nonlinear 
regime first. This is the basic picture of hierarchical structure formation: small scales 
collapse first to nonlinear structures, merging only later to larger structures. For 
WIMPs, this typically happens at a redshift between 40 and 80, depending on the 
interaction theory underlying the WIMP candidate and the details of the primordial 
power spectrum. All matter density perturbations that grow nonlinear are characterized 
by masses above Mf s . 

Whether these first purely gravitationally bound WIMP structures can survive until 
today is under debate. Their mean density contrast is comparable to that of the galactic 
disk and an order of magnitude above that of the halo in the solar neighborhood. One has 
to keep in mind, however, that the very first perturbations entering the nonlinear regime 
are, rather, rare fluctuations with mass scales M r f > Mf s . Those rare fluctuations can 
easily reach density contrasts which are an order of magnitude higher. High resolution 
simulations [26J show that structures with physical characteristics set by the cut-off in 
the linear power spectrum may survive merger processes of the first halos. The next 
critical phase starts with the formation of stars. Encounters with stars could disrupt 
the first CDM halos. So far, various studies [274 EE] show that the first halos would 
loose mass during an encounter, but their core should remain intact. 



5. Impacts for the indirect detection of dark matter 

Indirect detection methods for DM rely on the fact that most particle DM candidates 
can pair-annihilate into standard model particles. The annihilation products could then 
potentially be spotted in cosmic rays of various kinds, with an expected flux contribution 
that is proportional to the DM density squared. From hierarchical structure formation, 
as mentioned before, typical Milky- Way sized Halos are expected to contain a plethora 
of DM substructures rather than a homogeneous DM distribution. Focusing in the 
following on gamma rays, the consequences for indirect DM detection are two-fold and, 
as we will see, depend sensitively on an accurate determination of the cutoff scale in the 
power-spectrum as we have performed it in the preceding part of this work. 

At this point, let us briefly introduce some necessary notation and denote by 
p su b(M su b, r su b) the density profile inside a typical dark matter subhalo (or "clump") 
of mass M su b (as opposed to the smooth, i.e., average DM density p(r) in the Milky 
Way), and by n suh (M suh ,r) the number density of such clumps at a position r in the 
halo. For subhalos that cannot be resolved separately, the average annihilation flux 
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originating from r is then proportional to p 2 ff (r), rather than p 2 (r), where 

Pcfr( r ) = J dM suh J d 3 r suh n suh (M suh ,r - r sub )(p sub (M sub , r sub )) 2 (28) 

~ po J dM sub M suh S(M suh )n suh (M suh , r) . (29) 

The dimensionless quantity 

1 / d 3 r sub (p sub (M sub , r sub )) 2 
Po Jd 3 r sub p sub (M sub ,r sub ) 

that appears here is a measure of the effective density contrast between an average dark 
matter clump and the local halo density pq. The last step in ( |29l) is valid if n sub does 
not change very much on scales of the order of the size of a clump, i.e. for the whole 
region of integration one has n sub (M sub , r - r sub ) « n sub (M sub , r). 

Equipped with the above notations, we can now turn our discussion to the total 
gamma-ray flux from DM annihilations which is potentially strongly enhanced w.r.t. the 
case of a smooth halo distribution. To see this, note that the number density of subhalos 
scales with their mass as n sub oc M s ~ 2 - a behaviour that is expected from the theory 
of cosmological perturbations and has recently numerically been confirmed down to a 
mass scale of 4x 10 6 M Q , corresponding to the resolution limit of state-of-the-art iV-body 
simulations [Hj. To a first approximation, the density contrast 5, on the other hand, 
scales only weakly with M sub ; from ( |29i) . this means that the luminosity of individual 
subhalos is inversely proportional to M sub , leading to an equal contribution to the total 
subhalo-induced gamma-ray flux from each decade in subhalo mass - which is precisely 
the behaviour that has numerically been found in [H]. As we have stressed before, it 
is therefore crucial to have an exact and reliable method to determine the cutoff in the 
spectrum below which there are no contributions to the annihilation flux. Extrapolating 
the results of the above-mentioned study down to the typical cutoff mass scales that 
we found earlier, M& ~ 10~ 6 M Q , we expect a factor of almost 4 for the increase of 
the expected annihilation signal as compared to that from a smooth halo distribution. 
Taking into account the effects of sub-subhalos may further enhance the annihilation flux 
significantly; on the other hand, reduced survival probabilities for very small clumps, 
as discussed in the previous Section, would lead to slightly less optimistic predictions. 
Each of these issues is far from being settled at the moment of writing and warrants 
more dedicated future studies on their own. 

The other - perhaps even more important - effect of a large fraction of substructures 
in the DM halo distribution is the possibility of observing single DM clumps as gamma- 
ray point sources, maybe already in the near future with upcoming experiments like 
GLAST [29]. Comparing the projected all-sky GLAST sensitivity map for point-sources 
of DM annihilations [30] with the simulated gamma-ray sky of a Milky Way-like halo 
[14] . unfortunately seems to suggest that, in order to actually resolve the brightest 
substructures, point-source sensitivities about two to three orders of magnitude better 
than the GLAST performance are needed. One has to keep in mind, however, that 
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the flux estimates obtained in [2] are extremely conservative in the sense that no 
substructures smaller than ~ 10 4 M are resolved. This means, in particular, that the 
density profile within the clumps is essentially cut off at that scale; actual physical 
DM subhalos will exhibit much higher central densities and thus considerably enhanced 
fluxes (see also [Hj for similar conclusions). 

In order to quantify these last remarks, we adopt for the subhalo profiles the usual 
(a, P, 7) parametrization, 



and assume it to be valid down to a scale r cut = max{r kd , r ann } (here, r kd corresponds to 
the earlier determined cut in the power spectrum due to damping effects after the kinetic 
decoupling of WIMPs from the heat bath, and r ann is set by the maximal possible DM 
density due to DM self-annihilations [32]). The scale radius r s determines the transition 
from large (p su b oc r _/3 ) to small (p su b oc r~ 7 ) radial distances; for the small subhalos we 
are interested in here, one expects rather large concentration parameters [26l[3T], which 
means that the main contribution to the annihilation flux originates from the innermost 
regions of the subhalos and we can, for our purposes, well approximate the profile as 
p su b oc r~ 7 . Since the annihilation flux from a single subhalo is proportional to 5, as 
defined in Eq. (j3"0]) . we don't expect a strong dependence of this flux on the value of r cut 
for the case of mildly cuspy profiles (7 = 1). For a more cuspy profile like 7 = 1.5, on 
the other hand, we expect the annihilation flux to scale logarithmically with r cut , thus 
leading to an enhancement of up to 30 w.r.t. to the flux considered in [H]. In fact, 
at least at higher redshifts, one expects even steeper profiles for small subhalos, with 7 
ranging from 1.5 to 2 [26], thus potentially further enhancing the resulting flux. 

Of course, one has to keep in mind that neither the details of the innermost subhalo 
profiles nor the survival properties of these objects are sufficiently known at the present 
moment to make firm predictions. The numbers that we have presented above should 
therefore be considered as rather naive estimates, meant to illustrate the potential for 
GLAST - and certainly for successive, next-generation gamma-ray telescopes - to detect 
particularly bright subhalos as gamma-ray point sources. Note, in particular, that once 
GLAST has detected promising looking point sources, these could further be studied 
with air Cerenkov telescopes (ACTs) like CANGAROO [33J, HESS [33!, MAGIC [35] 
or VERITAS |36j . taking advantage of the much better sensitivity of these telescopes 
as compared to satellite based experiments (up to three orders of magnitude at high 
energies - see, e.g. [6]). 

6. Conclusions 

The dominant damping effect for the evolution of linear CDM perturbations is given 
by collision-less damping due to the geodesic motion of CDM particles. Even though 
collisional damping is sub-dominant, it has to be taken into account during the kinetic 




(31) 
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decoupling process, which provides consistent initial conditions for the free streaming 
regime. More importantly, collisional damping leads to a free streaming damping 
scale that depends on the temperature at kinetic decoupling. It is thus really this 
temperature scale, T kd , that renders the position of the cut-off scale in the linear power 
spectrum sensitive to the WIMP micro-physics. Since the cut-off scale might have 
observable consequences on sub-pc scales, at least in principal, the kinetic decoupling 
history provides an exciting interface between the WIMP interaction theory and the 
substructures that are probed by direct and indirect dark matter searches. 

We have presented, for the first time, a complete and analytic description of the 
kinetic decoupling process, without assuming a specific WIMP candidate; the master 
equation for the WIMP temperature evolution has been derived from first principles 
and accommodates any WIMP interaction theory. This could be achieved by separating 
slow and fast processes in the collision integral with a Born-Oppenheimer-type of ansatz, 
which proved to be an excellent approximation scheme for analyzing the WIMP radiation 
mixture - allowing, in particular, for a full error management; in this article, we have 
determined the decoupling scale to a precision of p 2 /M 2 ~ T kd /M x ~ 10~ 5 . 

The calculations presented here support the (order-of-magnitude) relaxation time 
approximation employed for the case of bino-like WIMPs before [9]; in this situation, we 
found a decoupling temperature that is about a factor of 2.7 smaller than the previous 
estimates, corresponding to a 4.2 times higher cutoff mass. As already mentioned, 
however, we were in particular able to significantly reduce the uncertainty in the 
determination of the cut-off scale. Let us stress again that this holds for arbitrary WIMP 
candidates and elastic scattering amplitudes. In particular, the method presented here 
can even be applied to regions of the parameter space where a simple scaling of the 
total cross-section with energy, e.g. o ~ uj° or o ~ uj 2 is no longer valid - as, e.g., in 
the case of s-channel resonances for LKP-like WIMPs. In this context, we would like 
to stress once more that it is, in fact, important to take into account the appropriate 
micro-physics; wrong assumptions about the particle nature of the DM candidate can 
lead to grossly wrong results in the determination of the cutoff scale. 

Turning to explicit models, we have performed a complete scan over the mUED 
parameter space and found a very narrow range for the kinetic decoupling temperature: 
27MeV < Tkd ^ 34MeV. A parameter scan over the MSSM with the methods 
presented here lies beyond the scope of this work, however, see [18] for a current 
state of the art approach (using the order-of-magnitude relaxation time approximation). 
Another interesting direction for future studies would be to develop a framework for a 
kinetic description even for DM candidates like axions that never have been in thermal 
equilibrium. 

Finally, we have commented on indirect detection prospects of a boosted signal from 
orthodox locations of enhanced dark matter densities, e.g. the galactic center or dwarf 
galaxies, and, much more promising, from individual dark matter clumps as gamma ray 
point sources. It remains, in fact, an exciting possibility that such objects may be seen 
already with the next generation of gamma-ray telescopes. 
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Appendix A. Boltzmann equation and temperature evolution in a 
Priedmann- Robertson- Walker spacetime 

The Boltzmann equation is one of the fundamental equations of non-equilibrium 
statistical mechanics and describes the evolution of the phase space distribution function 
/(x, p, t) for a given particle species \- m an abstract form, it is usually written as 

L[f) = C[f) . (A.l) 

Here, the Liouville operator L is the covariant generalization of the convective derivative 
familiar from hydrodynamics, or - in more technical terms - the variation with respect 
to an affine parameter A along a geodesic: 

L\f\=# = + = rf^- - I* (A 2) 

lJi dX dX dx» dX dp» y dx» paP y 0p» ' 1 ' ; 

Since x is supposed to be on shell, p° is actually not an independent degree of freedom - 
so in the above equation, we have df /dp = 0. Let us now consider a (flat) Friedmann- 
Robertson- Walker (FRW) spacetime, 

ds 2 = dt 2 -a 2 (t)dx. 2 , (A.3) 

and express everything in local (or comoving) momenta p = ap: 

= P° U + a- 1 p ■ V + ^p ■ V p -ar^yAj / (A.5) 

= p°(d t -Hp-V p )f. (A.6) 

In the last step, we have introduced H = a/a and used the fact that the homogeneity 
of space implies / = /(|p| , t). From now on, we will only use comoving momenta; for 
convenience, however, we will in the following drop the bars over the p. 

The collision term C[f] on the right hand side of the Boltzmann equation takes 



account of all interactions that x may experience. In Appendix B we determine its 
exact analytical form for a non-relativistic particle that is in contact with a heat-bath 
of relativistic particles F through two-body scattering processes of the type Fx — > Fx- 
As explained in Section [2j this is precisely the situation that we encounter after the 
chemical decoupling of x i n the early universe. To leading order in T/M x and p 2 /M 2 , 
where T is the temperature of the heat bath, the collision term then reads 

(rp \ n+4 
—J {M x TA p + Ap.V p + B}f(p). (A.7) 
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In the above expression, n is determined by the low-energy behaviour of the scattering 
amplitude, see Eq. flB.211) ; c is an overall constant that is given in Eq. flB.22j) . 
Furthermore, we actually have A = 1 and B = 3, but in order to develop a physical 
understanding for the value of these constants, we will, for the time being, keep our 
discussion to the general form displayed in Eq. (1A.7I) . 

Let us start by dividing the Boltzmann equation by p° and then integrate it over 
all momenta. On the left-hand side, we find 

d 3 p 



f27r) 3 



R-#p-V p }/(p) 



< rn x) 



(A. 



where n x is the number density of \. Performing the same integral over (1A.7|) gives 

73^ / m \ "+4 



drp 



C[f]/p° = cM x 



T 

mZ 



(B - 3A)n 



x ■ 



(A.9) 



Of course, scattering processes should not change the comoving number density a 3 n x - 
which, from the above expressions, is only true for B = 3A. 

Inspired by the situation for a thermal distribution, let us now define the 
"temperature" T x for an arbitrary distribution / of non-relativistic particles as 

d 3 p 



P 2 /(P) 



3M X T x n x . 



(A.10) 



With this definition, we have T x w 



T for a situation close to thermal equilibrium and 
e = (T — T x )/T can therefore be considered as a parameter that characterizes the 
deviation of / from a thermal distribution. To determine the evolution of T x , we consider 
the next moment of the Boltzmann equation: 



(2vr) 3 p° 



(A.ll) 



Performing the integrals, and using B = 3A, we arrive at the following process equation: 



with 



T' - 

x 



a = 



n+2 



T 



X 

T 



T 

Ml 



n+2 



45 

7T 3 ^ e fr 



1/2 



Mi 



pi 



(A.12) 



(A.13) 



Here, T' x = dT x /dT and we have assumed that the universe is radiation dominated, with 
g e g effective degrees of freedom in relativistic particles. For n ^ —2, Eq. (IA.12I) can be 
solved by the ansatz 



T x = A(T) T x hom) 



A(T) 

M v 



T 2 exp 



aA 



n 



T 



2\M V 



n+2 



(A.14) 



The general solution to the inhomogeneous differential equation (1A.12j) then reads: 

z l/(n+2) 



T 

T = — 

x A 



1 



n 



exp[z]r[-(n + 2)- 1 ,,]) N , l+2 + A T( hom ). 

) z— aA I -2L I 



(A.15) 



Thermal decoupling of WIMPs from first principles 



21 



535 
o 

bJ3 
O 




logio x 



Figure Al. Phascplot for the temperature evolution, expressed in the dimensionless 
quantities x = \a\ (T/M x ) n+2 and y = T x /T. Here, we show the case for n = 0, but 
the same type of behaviour is obtained for any other (positive) value of n. As can be 
seen, any departure from thermal equilibrium (T x = T) is restored almost immediately 
for x > 10, forcing T x to follow the solid line which represents our general solution 
(|A.22[) . For smaller temperatures, \ decouples from the heat bath; in this regime, its 
temperature is only influenced by the Hubble expansion and thus evolves as T x oc T 2 . 



In the limit of large and small T, respectively, the incomplete Gamma function behaves 
as 



(n + 2 



-1 



z] T -=°? z-V(«+2) exp 



(n + 2) 



n 



r[- 

r [" V'" ' "J J ~J ' ^l/(n+2) 

Only for A = 1, as obtained in (IB. 221) . one can thus recover the expected asymptotic 
behaviour (by choosing Aq = 0): 



z]/z, 
r[-l/(n + 2)] . 



(A.16) 
(A.17) 



T X T -^T 



n + 2 



l/(n+2) 



n + 1 
n + 2 



(A.18) 
(A.19) 



At high temperatures, the particles are still kept in thermodynamic equilibrium, 
i.e. T x = T. The scaling T x oc T 2 at low temperatures, on the other hand, corresponds 
to the solution of the collision-less Boltzmann equation, i.e. Eq. (1A.12|) for a = 0. 
To conclude our discussion of the process equation (1A.12[) . let us rewrite it as 



(n + 2)^ = y±(y-l), 
ax x 

where we have introduced dimensionless variables 



x = \a\ [T/M x ) n+2 



and 



v = t x /t, 



(A.20) 



(A.21) 



and the upper (lower) sign applies if a > (a < 0). In this form, we can immediately see 
that the asymptotic solution y — 1 for x <S 1 (which corresponds to the expected T x = T 
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for large T) is stable against small perturbations if and only if a > 0, i.e. c > 0. In fact, 
from the phaseplot shown in Fig. IA1I it becomes apparent that in this case \ is kept 
in thermal equilibrium extremely efficiently for x > 10 and that therefore the evolution 
of T x actually follows an attractor solution to Eq. (1A.12I) . Once the temperature falls 
below x ~ 10, decoupling occurs on a rather short timescale and we observe the expected 
T x qc T 2 due to the Hubble expansion. 

To summarize the results of this Appendix, we have developed a physical 
understanding for the functional form of the collision term ( 1A.7I) and found that T x 
evolves as a function of T according to 

( -,1/0+2) } 

T x = Tll-— TY eMz}T[-{n + 2)-\z\\ , (A.22) 

2— n + 2 \ M x ) 

essentially independent of the initial conditions. For convenience, we state here again 
the full expression for the constant a that appears above: 

/ 10 \ M 

a = 9suc n Nt +3 — , (A.23) 

where the sum runs over all relativistic scattering partners / (with g^M internal degrees of 
freedom) and c n is the corresponding scattering amplitude at zero momentum transfer, 
see Eq. ( IB. 211) ; g e g is the effective number of heat-bath degrees of freedom close to 
decoupling and N~f is defined in Eqs. ( IB.23[ IB.24I) . Finally, we may now make use of 
the asymptotic behaviour, Eqs. flA.18llA.19l) . of our general solution to define the kinetic 
decoupling temperature by equating these two limits: 

-i 



rp / rp 2 \ I / n \ l/(n+2) 



M x \M X TJ T ^ \\n + 2 



r 



n + l 
n + 2 



(A.24) 



Appendix B. The collision term for two-body scattering processes 



In this Appendix, we provide the details for the calculation of the collision term for 
scattering processes between a non-relativistic particle \ an d relativistic (usually SM) 
heat-bath particles. For the former, we use p^ = (E, p) to denote ingoing momenta, 
while for the latter we use = (u, k); the corresponding outgoing quantities are marked 
with a tilde. With these notations, the collision term takes the form 

where Ai is the scattering amplitude, summed over final and averaged over initial spin 

states, and 

J = 9SM [(l T 9^(00) g+mip) "(IT g^) g±{v)f(p)] ■ (B.2) 

We will assume that the relativistic SM particles are thermally distributed, with (?sm 
being the number of spin states, so their phase-space number density is given by 



^(^±1)^ (B.3) 
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for fermions (upper signs) and bosons (lower signs), respectively. The only assumption 
that we will make about the \ distribution function /(p) is that we can neglect Pauli 
suppression factors - as we have done in ( IB. 21) . In the context that we are considering 
here, this is satisfied to a very good approximation. 

Since, for kinematic reasons, the average momentum transferred during the 
scattering events is small, we can expand fIB.lj) as 



C{E)= 



(B.4) 



d 3 k 



(2tt) 3 2u;J (2vr) 3 2cD J (2vr) 3 2£ 



d 3 p 



1 



J! 



^(P)5 (3) (P"P) 



(B.5) 
(B.6) 



x (2n) A 5(E + Cj - E -u)\M\ 2 J 

where we have introduced 

D q (p)=q-V p = (k-k)-V p , 

and the derivatives of the delta function that appear above are as usual defined in terms 
of integration by parts. 

After these preliminaries, let us now start to calculate the expansion coefficients 
. First, we note that J\p =p q =u = 0. From this, it immediately follows that 

C° = 0. (B.7) 

For the next coefficient, we find 
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-2tt 
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d 3 k 



(2tt) 3 2cj J (2it) 3 2Cj 
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\M\ 2 J 
2E 



S(E + Co - E - u) 



(B.8) 



p=p 



7T 



d 3 k 



d 3 k 



x 5 (to — to) 
where we have used 
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(B.9) 



(B.10) 



to perform the necessary integration by parts in Co. Let us now write the scattering 
amplitude as a function of the Mandelstam variables s and t, where s = M 2 + 2u(E — 
|p| cos9), and expand it in |p| / E: 

^2 



\M 



t=o 



\M\U 

e=Tr/2 



uj-^cos9 [ d u \M\ t 



2 

t= 

0=71-/2, 



(B.ll) 



With 1 || p || V p J\p =p , the angular integrals that appear in the above expression for 
C 1 are then of the form 



16 



dn j dn \M\t =0 (c L -\) = f7T-WE- 1 (p-\){ \du\M\l, 
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To leading order in p 2 /E 2 , and with the help of the useful relation 



g ± (lT9 ± )(uj) = . = 

(e«/T ±1 )2 

we thus arrive at the following result for C 1 : 

fi'SM 



(B.13) 



Cl = -^ M x 3 J duu^d^) U\M\lA [™ xP -V p + p 2 ] /(p). (B.14) 



Let us now compute the next coefficient of the expansion (1B.4I) : 
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p=p 



This time, the angular integrals are of the form 

dQ J dQ \M\ 2 t=0 gW = y^ 2 (u> 2 + 2 ) \M\l y j +O (j^j .(B.17) 

Again, we will only keep the leading order terms and, after similar steps as for the 
computation of C , we arrive at: 



C 2 



fi'SM 

6(2vr) 3 

x 



M- 3 f duu* \M\ 2 t=0 

J e=n/2 



[-TM 2 d UJ g ± A p + M x (ATu^d^ - d u g± + Td 2 w g±) p ■ V p 
-SM x d u g± + p 2 (Auj-^g* + d'g*)} /(p) . (B.18) 

The terms C- 7 for j > 2 only contribute higher-order corrections in p 2 /E 2 and 
uo/M x . To lowest non- vanishing order, and after the necessary partial integrations, we 
thus find the following expression for the collision integral: 

C = C l + C 2 (B.19) 

fi'SM 



6(2tt 
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(B.20) 



Unless one encounters s-channel resonances (see Appendix C for a treatment of this 
case), one may now expand the amplitude at zero momentum transfer in uj/M x : 

n I / \ n+l N 



I A* 



i=0 



+ 



(B.21) 
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Figure CI. Feynman diagrams for the elastic scattering of the LKP with SM fermions. 
In principle, there is also a contribution from a i-channel Higgs-exchange, but this is 
strongly suppressed due to the small Yukawa couplings involved. 



Keeping only the lowest order term of this expansion, even the integral over the photon 
energies can be performed analytically and we arrive, finally, at the following expression: 



where (for j > 0) 



.7 + 1 
.7 + 1 

rpj + 1 



"n+3 



M X TA P + p 
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(1-2-^ (j + l)!C0' + l), 
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Hp), (b.22) 

(B.23) 
(B.24) 



Note that in (IB. 221) - as in all corresponding formulas before - a summation over all 
SM scattering partners is understood. 

The expression (1B.22j) is the central result of this Appendix; let us conclude by 
remarking that the formalism which we have presented here in some detail can straight- 
forwardly be used to calculate even higher-order corrections to C. However, as we have 
learned in Appendix A x wm t> e kept in thermal equilibrium until p 2 /E 2 , (u)/M x ~ 
10~ 6 — 10 -5 , so the approximation to keep only leading order expressions in these 
quantities guarantees a result that is correct to a very high accuracy. 



Appendix C. The LKP scattering amplitude 

In this Appendix, we present the LKP scattering amplitude close to kinetic decoupling, 
which is dominated by the scattering with light SM fermions. The relevant Feynman 
diagrams are shown in Fig. ICll Note that in the UED model, at each given KK level, 
the number of fermions is doubled as compared to the SM, i.e. singlet and doublet states 
come each equipped with a whole KK tower of their own. The corresponding couplings 
of the LKP, well approximated by the are given by (see, e.g., [37]): 

- VJV(1 + 7 5 )/ (0) + cc. (C.l) 
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0.0001 0.001 0.01 0.1 1 



u/M B (i) 

Figure C2. The full LKP scattering amplitude (solid lines) vs. the approximative 
expression that scales with the photon energy as u> 2 (dashed lines). For simplicity, 
we have shown the case of electron scattering, and only taken into account the 
contribution from left-handed states. The two situations correspond to a mass splitting 




/M S (i) of 1% (left) and 5% (right), respectively. 



for KK fermion singlet states, /i , and 

for KK fermion doublet states, f^\ Given these couplings, and after summing over 
all initial states, it is straight-forward to calculate the amplitude at zero momentum 
transfer (in the limit of rrif=0): 

2 



(C.2) 



\M\l 



s=M%+2M x u 



+ i6M2 (1) M; (1) ^ 2 r 2 

Is 



+ (s^d) , (C.3) 

where T is the decay width of the corresponding s-channel state. For fermion energies 
much smaller than the mass difference between the LKP and the KK leptons, w/M 5 (i) <C 
5 = [Mf(i) — M B (i)j /M B (i), we may expand the scattering amplitude as we did in 
(E2J]): 
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In Fig. IC21 we compare the full expression for the scattering amplitude with the 
approximation of keeping only the lowest order term in u/M B (x>\ we can see that the 
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low-energy approximation is valid as long as u < 0.15M B (i), and then starts to deviate 
considerably. Considering that we expect decoupling temperatures of ~ 10 MeV (see 
Section ISTBl . this condition is easily satisfied for the average photon energy, (u) = 2.70T, 
during the decoupling process. Still, one might worry about photons in the tail of the 
thermal distribution, which may very well hit the resonance and thereby enhance the 
effective scattering amplitude considerably. 

To address this question, we have to go back to the full expression flB.201) for the 
collision integral that we derived in Appendix B Going through similar steps as those 
leading to (1A.12|) . we then find that the process equation for T x to be solved is 

r x -[2 + a(T)]T x /T=-a{T), (C.6) 



where 



< T ) = E (to^t)^ ^mJ^ 2 f dug+WiK, L 4 |ACo ) (C7) 

and the sum runs over all relativistic scattering partners. We may now compare a(T), 
for T ~ T kd , between the two situations where we either take the full expression (1C.3I) 
or the approximate expression (1C.5I) for the scattering amplitude. Using 

T = ^6 2 M BW , (C.8) 

which is valid in the limit of small mass splittings 5, we then find that for 

5 > 3 x 1(T 3 (C.9) 

the mistake in a(T) is less than 5% - which, from (IA.24I) . roughly, corresponds to a 1% 
error in the estimate for T^- For smaller 5, on the other hand, the two ways to calculate 
a(t) soon start to give grossly different results, rendering this type of error estimation 
invalid; in such a situation, one has to (numerically) study the asymptotic solution to 
the full process equation (1C.6I) in order to determine the decoupling scale. This shows 
that, indeed, the high-energy tail of the heat-bath fermions plays an important role in 
such a situation. 

Recalling that in the mUED model we have 5 ~ 10~ 2 , we will here not further 
analyze the case of a very pronounced s-channel resonance. However, we note that 
mass splittings of the order of 5 ~ 1CT 3 (and smaller) are actually feasible in, e.g., 
non-minimal UED versions. Here, we have provided all the necessary tools to correctly 
treat such situations. 
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